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Abstract 

The theory describing action of medicines explored in this paper is based on as- 
' ^ . sumption that vital activity of the cell may be described in terms of the model of two 

r1 ■ states: resting state and excitation. According to available physiological data excita- 

tion state is dangerous for cell and may cause different pathological changes, including 
"conformational" diseases, due to protein aggregation. Normally, the excitation is com- 
pletely reversible and the key role is played here by ATP (adenosine-5 '-triphosphate) 
which disaggregates proteins of cytomatrix. The same effect ATP exerts during cell 
injury by eliciting a "healing" effect. Damage of cell structures we consider as "ill- 
ness", whereas removal of pathological consequences caused by protein aggregation of 
any origin we will call "a cure". The latter is considered as physical process of cell 
recovering from excitation/injury to resting state, which is analyzed in terms of our 
generalized thermodynamics. "Cure" results in reduction of effective temperature of 
cell proteins due to binding intracellular water which is present in the cell at concen- 
tration of approximately 44 M. As a result of the sorption processes mobility of both 
water and proteins is greatly reduced, and the corresponding reduction of the effective 
temperature seen by us as a measure of treatment effect. 
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1 Introduction. 



The data collected by classic cell physiology revealed that delayed excitation phase leads to 
damage of cell structures (Nasonov, 1962), followed by different cellular pathologies, which 
we call "illness". In order to perform physical analysis of "illness" and "cure" one should 
rely on an effective model of cell which, in this case, we have chosen to be Ling's model 
(Ling, 2007). The principle parameters of this model are based on assumption that vital 
functions of cell are described as cyclic transitions between two states: resting state and 
excitation (two state model). According to Ling, a key factor which "switches" cell between 
these two states is ATP. At resting state ATP is adsorbed on key cell proteins, that leads 
to formation of the following complex: (ATP)m(protein)„(H20)p(K"'")q, where in may be so 
big so that all key cell proteins may become a part of such complex; complex with n = 1 
is considered as physiological atom or physioatom (Matveev, 2005), as in this case complex 
under consideration represents a minimum structure which preserves capacities of a whole 
resting cell. 

Physioatom at resting state is characterized by condition when protein inside the complex 
has an unfolded conformation so that all its peptide groups are available to the solvent and 
can absorb water creating from water molecules a multilayer relatively ordered structure 
(adsorption layer), whereas free carboxyl protein groups selectively bind K+ ions in the 
presence of Na+ (Ling, 2007). 

During ATP enzymatic degradation physioatom is being destructed leading to appearance 
of the following products: m(ADP) + n(Pj) + p(H20) + q(K^), where Pj - orthophosphate. 
In other words, water and K+ are being desorbed, whereas protein molecule folds so that its 
peptide bonds become unavailable to the solvent. At this state proteins have high aggregation 
activity, which is responsible for appearance of novel protein structures (Matveev, 2010). 
Altogether, these changes transit cell to excitation state. Inflow of new ATP molecules 
inhibits excitation process and eventually returns cell back to the original state - resting 
state. 

It is well known that delayed excitation phase cause damage of cell structures that is 
a reversible process at early injury stage. While irreversible damages keep increasing cell 
dies (Nasonov, 1962). For short, reversible damage phase will be called as "illness", where 
influences that contribute to cell recovery (return back to resting state) will be known as 
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"healing", and process of cell recovery - as "cure". 

According to Ling's model of the cell the most effective "medicine" is ATP, however 
other compounds may have the similar, in principle, effect on proteins; such compounds we 
will also refer to drugs. Here we present physical analysis on mechanism of action for such 
medicines. We will make use of so called generalized thermodynamics which we proposed 
earlier (Prokhorenko and Matveev, 2011). By using theoretical analysis we were able to 
obtain the following results, (i) Transition of the cell to excitation state is an exothermic 
process, i.e. heat is being discharged under destruction of physioatom. (ii) During transition 
from resting to excitation state cell volume is being reduced, (iii) There has been proved 
principle of Ling's theory stating that key cell proteins undergo folding under excitation, 
(iv) By using analytical approach there was shown that K+ ions are indeed being freed from 
bound state followed by release from excited cell along concentration gradient. 

By using approaches which we proposed there was possible to set up a task of formu- 
lating physical theory for mode of action of medicines resembling ATP (Ling's medicines). 
Moreover, it is obvious that any influence on damaged cell will be considered as "healing" 
effect if they are able to increase intracellular ATP content or lead to appearance inside cell 
of any other compounds with ATP-like properties. On the other hand, our description of cell 
injury may be considered as physical theory of cytopathology, and more special — physical 
theory of acting poisonous compounds. 

According to Ling, effect of ATP on sorption properties of protein is explained by its 
electron acceptor capacity: when ATP binds to protein it causes decrease of electron density 
on functional groups of protein so that it affects their possibility to interact with other 
functional groups of protein, water, and major cell cations K+ and Na"*". Such shifts in 
electron density result in destruction of secondary protein structures as well as to increased 
affinity of peptide backbone for water; free carboxyl groups of protein will begin to bind K+ 
in the presence of Na"*". 

Such statements of Ling's theory we consider as principles so that their physical nature 
will not be a goal of the current study. Instead, we will point out at chemical compounds 
(medicines) that act on the relevant proteins mainly resembling that one for ATP. If chemical 
modification of protein leads to the same changes as if they were induced by ATP adsorption, 
such influences will also be considered as healing. As a result of analysis we were able to 
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come to the idea about effective temperature of physioatom as a measure of healing effect: 
while physioatom is being developed substantial amount of water molecules were bound, 
then motility of both water and physioatom as a whole will be markedly decreased. Overall, 
stabihty of physioatoms themselves as well as their associates will be growing up, thus 
enhancing their resistance to perturbing factors of external and internal milieus including 
thermal influences. 

The current paper is structured in the following way: part 2 contains overview of gen- 
eralized thermodynamics; part 3 is dedicated to description of physical theory for mode of 
action of some medicines; part 4 is aimed at discussion of ATP as a Ling's medicine in the 
context of the paper; part 5 - conclusion. 

2 Overview of generalized thermodynamics 

Absence of models for living cell which might be compatible with the routine approaches of 
theoretical and mathematic physics represent a substantial obstruction for introduction of 
physical methods into biology rather than an illusive gap between animate and inanimate 
natures. To the best of our knowledge, only one such model is in the literature - a model 
of hving cell proposed by Gilbert Ling (2001), which we investigated from the point of 
thermodynamics and statistical physics (Prokhorenko and Matveev, 2011a,b). 

We showed that after performing some natural generalization conventional equilibrium 
statistical mechanics then may be successfully applied to explain a number of events which 
take place in Ling's cell when it gets activated or dies. Such generahzation is based on 
the property of nonergodicity, which was recently proved for Bose gas with weak pairwise 
interactions (Prokhorenko, 2009). However this proof was so general that it was true not 
only for Bose gas but for the majority of the most realistic systems of statistical mechanics 
as well. Property of nonergodicity states that apart from trivial motion integrals (first 
integrals), e.g. Hamiltonian (energy), impulse, number of particles, there are the others 
commuting independent first integrals, which commutes with trivial ones. 

Approach that we suggested within equilibrium statistical mechanics and thermodynam- 
ics is called generalized thermodynamics (Prokhorenko and Matveev, 2011a) which is based 
on primary elements such as: (1) a number of motion integrals of system commuting with 
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each other (in involution), (2) generahzation of the standard Gibbs microcanonical distri- 
bution for occasions when the first integrals from the mentioned commuting integrals take 
up defined values, and (3) entropy as a function of values for energy of system and first 
integrals. 

Let us consider the main elements of generalized thermodynamics in more detail. Let 
us assume that apart from Hamiltonian H there are also / = 1,2,3... integrals Ki,...,Ki 
which commute with each other. Property of integrals Ki,...,Ki to commute with each 
other is necessary to simultaneously measure them, thus, attributing to them defined values. 
Detailed analysis of measurement was performed by von Neumann (1932). 

The proposed generalized thermodynamics is built upon the equation for distribution 
function in phase space (statistical matrix) given that energy of system takes up value E 
and integrals Ki, ...,Ki - values K[, respectively: 



Entropy of system that corresponds to this distribution (in case of classic mechanics) is 
shown as: 



where dVx - is an element of phase space. In quantum mechanics integrating over phase 
space shown in the latter equation must be replaced by operation of taking trace. Equations 
([T]) and ([2]) represent generalizations for microcanonical distribution and its entropy in case, 
when system contains non-trivial (commuting with each other) first integrals. 

Let us answer a question: when do equations ([1]) and ([2]) may give novel results if they 
are used instead of microcanonical distribution? For simplicity let us consider an occasion 
when 1 = 1. Then, under specified energy of system E entropy S{E) is a function of value K 
of first integral K' . There may be two most typical dependencies for value of this function 
on K': 

1) S{E, K') has maximum on K at isolated point K' = K. In this case equilibrium will 
take up a value K' = K so that, as has been shown before (Prokhorenko and Matveev, 
2011a), it will be possible instead of distribution (1) to use a regular microcanonical distri- 
bution. In other words in this case if one uses generalized microcanonical distribution no 
benefit will be obtained compared to using standard microcanonical distribution. 




(1) 




(2) 
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2) S{E, K) reaches maximum at some interval of non-zero length, i.e. S{E, /T) as a func- 
tion has a plateau. In this case generalized microcanonical distribution will give results that 
will be different to those obtained after applying standard microcanonical distribution. This 
is not something special. In theory of phase transitions it is quite often when thermodynamic 
values have plateaus (as function of parameters of system). The main assumption is that 
(Prokhorenko and Matveev, 2011a,b) the latter of two occasions corresponds to Ling's cell 
at resting state. 

Generalized thermodynamics has many similarities with the standard thermodynamics. 
For example, the standard formula dE = TdS — PdV is also true for it (T is temperature, 
P is pressure, V volume of system). However, for generalized thermodynamics some other 
new results are true as well, which is due to the presence in system of these non-trivial first 
integrals commuting with each other. Let us describe them. 

Commuting first integrals Ki, ...,Ki which are involved in formula for generalized mi- 
crocanonical distribution (1) are called active. Note that, not all integrals among some 
maximum set of independent commuting integrals but only active ones may be used in 
distribution (1). 

We think that for Ling' cell there are so many independent commuting first integrals 
so that a number of active integrals may be described by some continuous parameter s e 
[0, 1], given that s = corresponds to the absence of active first integrals, whereas s = 1 
corresponds to active integrals, which constitute a maximal set of commuting independent 
first integrals. Existence of continuous parameter s, that "numerates" first integrals of 
system commuting with each other has been explained earlier (Prokhorenko and Matveev, 
2011 a. Let us start to describe main results of generalized thermodynamics. 

I. let, s ^ s~ ds, 6s - is an infinitely small quantity. Let {SS)e - is entropy change under 
this process when energy does not change. Then, {SS) e is constant under adiabatic change 
of system variables i.e. 



This statement has been proven (Prokhorenko and Matveev, 2011a). 
II. Again, let us assume that s ^ s — 5s, 5s is an infinitely small quantity and {5S)t-i 
{5S)e), is a change of entropy and energy of system in this process under stationary tem- 




(3) 
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perature. Then 



:(55)t<0. 



This statement was proven as well (Prokhorenko and Matveev, 2011a). 

III. Finally, to explain physiological events in Ling's cell additionally we used another 
obvious principle: there is such temperature above which life can not exist. 

Above such temperature entropy as a function of the first integrals Ki,...,Ki will loose 
plate corresponding to maximum entropy, so that generalized microcanonical distribution 
([1]) reduces to the usual equilibrium distribution. In our previous work (Prokhorenko and 
Matveev, 2011a) this statement was not proven, but rather illustrated by several analogies 
taken from theory of phase transitions. 

Described above main propositions of generalized thermodynamics let us to explain 
the following physiological events which take place under activation and death of the cell 
(Prokhorenko and Matveev, 2011a): 

1. Excitation and cell death are exothermic processes. 

2. Size (volume) of cell is not kept stationary (mainly it decreases) when it gets activated 
or dies. 

3. When Ling's cell gets activated or dies key cell proteins are being folded. 

4. Activation or cell death are accompanied with efflux of K"*" ions out of cell. 

In the latter case we were able to link efflux of ions with cell heat release under 
excitation or death of cell. Recently we proposed two microscopic models of protoplasm 
(Prokhorenko and Matveev, 2011b) that might allow us to test statements of generalized 
thermodynamics in some particular numeric examples (Prokhorenko and Matveev, 2011a). 
The first one we call van der Waals model, which is characterized by the following features: 
1. At resting state all proteins are located at the sites of some crystal lattice. 2. It is 
postulated that interaction forces among proteins (van der Waals forces) are small enough 
to be considered protoplasm proteins as an ideal gas. 

Basing on data of cell heat release that take place under its excitation and death we 
showed that when cell dies (under transition to equilibrium state) all proteins which were 
unfolded now start fold and aggregate by forming 220 superclusters, which is controlled 
by non-covalent short-range forces (compared to van der Waals forces), but such forces are 
great in magnitude (physical and biological meaning of superclusters that were obtained 



theoretically are unknown yet). According to the proposed equilibrium conditions of protein 
aggregates (Prokhorenko and Matveev, 2011b) we calculated that the volume which proteins 
occupy in dead protoplasm is 0.136 (volume of dead cell is taken as 1), that correlates well 
with empirical value 0.163 calculated (Prokhorenko and Matveev, 2011b) on the basis of data 
regarding (mass) content of hemoglobin in (per volume unit) (Van Bcckvclt ct al., 2001) as 
well as data on molecular volume for folded hemoglobin (Arosio et al., 2002). 

In our second model we consider protoplasm at absolute zero temperature and a priori 
assume that all protein molecules are located at the sites of some crystal lattice keeping into 
consideration inner degrees of freedom for protein molecules and specially emphasizing on 
study of conformational changes in proteins. We have shown that the ground state of the 
second model (within the framework of some approximation) is the ground state of some 
Hamiltonian, which describes superfiuid Bose gas on some complex configuration space of 
protein molecule rather than in Euclidean 3D-space (as it is in the first model). By using 
this effective Hamiltonian and our understanding of living matter we were able to show that 
when Ling's cell is at the living resting state some of its proteins are unfolded whereas in 
the dead protoplasm all proteins are folded. 

Additionally, we consider a special occasion when cell transits from living state to equi- 
librium or dead (death process) state that happens in some weak external field, e.g. elec- 
tromagnetic or field of acoustic oscillations (which is a sum of harmonic motions with a 
continuous frequency spectrum and random independent phases evenly distributed over cir- 
cumference). For this process we derived the Fokker-Planck type equations, which describe 
spreading of distribution function at (reduced) phase-space of cell, and showed that external 
field in this equation comes only trough spectral density of intensity for zero frequency. 

This analysis is of importance because it allows to experimentally test the correctness 
of our calculations because the spectral field intensity corresponding to frequencies differing 
from zero is unable to affect state of a cell, as it follows from our conclusions. 

Compatibility of Ling's cell with our analytical approach admit us to investigate an issue 
of injury in Ling's cell as well its recovery (healing) after applying pharmaceutical agents of 
a new kind, which we call Ling's medicines. 
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3 Physical theory for mode of action of medicines 



In the current part we will present our physical theory of action of Ling pharmaceutical agents 
on living cell assuming that Ling's model of the cell is amply realistic. For simplicity, here 
and after, we will discuss our ideas in a classic descriptive way (i.e. within classic mechanics). 
We suppose that the living cell 9^ represents a Hamiltonian system (from classic mechanics), 
which is described by Hamiltonian H{p, q), where p and q are canonically conjugate variables 
of phase space of a cell, p are canonical momenta, q are canonical coordinates of cell. However, 
for us such description of the cell is supposed to be rather generalized, so we would like to 
consider dynamics of protein molecules separately. It turns out to be possible due to so 
called adiabatic limit, which we can use because of the very high mass of protein molecules. 

Thus, let us assume that cell volume is V and V is a volume for some selected area 
D, which is defined in way so that all flows of all particle through the boundary is equal 
to zero. Let Pq, q^ be canonically conjugate coordinates and momenta for protein molecules 
inside domain D which describe inner degrees of freedom for protein molecules. Let p'i,q'i, 
be canonically conjugated coordinates and momenta for the rest of particles present inside 
area. Thus, p' = {p'q,p'i) and q' = {q'Q,q'i) describe canonically conjugate coordinates and 
momenta inside area D. 

Likewise, let pg are canonically conjugated coordinates and momenta, of particles 

in addition to area D (for shortest it will be called as ). Let p'( and q'( be canonically 

conjugate coordinates and momenta for other particles from this addition to area D. Thus 
p" = {PqjP'D and q" = {qQ,q'l) describe canonically conjugate coordinates and momenta of 
particles from '^K\D. Variables Pq, satisfy to canonic Hamilton's equations 



Let us assume that cell has temperature T. As motion of proteins is very slow we can 
average right part of equations (jlj) on Gibbs conditional distribution for the whole molecules 
in domain D provided that momenta and coordinates for protein molecule in area D are 
equal to p'^ and q'^. Application of such average is called an adiabatic limit. After applying 



dH{p, q) 
dq'o 
dH{p, q) 

dp'o 



(4) 
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this average an equation will be as: 



ry, , /X / dp^dq^dp dq —— — e t 



and by definition 



Z{Po,qo) ^ J dp[dq[dp"dq"e t'"^ 



We will consider the domain D, although small, however, being a macroscopic so that in- 
teraction between particles in domain D and outside it may be neglected. In other words, 
complete Hamiltonian H(p, q) for cell may be considered equal to a sum of some Hamiltonian 
H'{p', q'), which describes a part of cell inside the domain D, plus Hamiltonian H"{p", q") , 
which describes a part of cell outside area D. In other words, one may conclude that 



H{p,q) = H'{p\q') + H"{p\q") 



Let us calculate 



1 /",',' dH'ip', q') H(r>',,') 

where by definition 

Z'iPo,%) ■■= J dp\dq[e 

Assuming that 

i^D(p'o,g^) = -TlnZ'(p[„g^). 

One can check that in adiabatic limit evolution of system composed of protein molecules is a 
Hamiltonian evolution corresponding to Hamiltonian Foip'o, Qo)- other words, in adiabatic 
limit there will be: 

dFoip'o^q'o) 



dq'o 



~ dp', 
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Of note, due to canonical Gibbs distribution for the whole cell density function for momenta 
and coordinates of protein molecules in area D can be described by distribution function 

p(Po,%) = ^e"^^T^, (5) 

where Zjj is a normalization factor 



Zd= j dp^dq^e t 
To simplify further calculations put by definition 

Now it is correct to ask a show does it happen that temperature T of protein molecules T', 
which is a part of distribution function for proteins ([5]), is equal to temperature of the whole 
cell T. The evolution of protein system is a Hamiltonian (in adiabatic limit) and can be 
described by Hamiltonian G(pQ,gQ), thus, being a motion integral. Thus, it seems evident 
that it should be T' = T. T' might be an arbitrary prescribed value corresponding to any 
given in advance value G' for integral G{p'q, Qq) so that the following equation will be 

dG' ' 

where 



S{G') = \nj dp',dq',6{G{p',,q'o) - G'). 



Of note, however, that we did not take into account here the following two issues: 

1. Fluctuations of parameters of domain D (e.g., its volume). 

2. Energy dissipation of protein molecules in area D. 

Accordingly, Liouville equation for distribution function p of protein molecules in area D 
may be calculated as: 

d 

—Pt = Lfpt + Ldpt 

Here, is a part of complete Liouville operator for protein system in domain D, that puts 
into account evolution of system due the Hamiltonian G{p'q, q'^) (corresponding to fluctu- 
ations of parameters of domain D), and is a part of Liouville operator of our system. 
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corresponding to dissipations that was unaccounted in G(pQ,gQ). Firstly, let us consider 
equation 

|a = LfPt (6) 

For this, we will assume that G{p'q, q^) only depends on the volume of system V, and instead 
of G we will write Gy- Let Vq be a time average of volume V, V Vq identically as the 
volume of area D may fluctuate. We suppose that volume fluctuations . are small, 

Gy may be expanded in a Taylor's series by V and we can rest only first-order terms. In 
other words, one may take 

G = H + 6{V-Vo)P, 

where H = Gvo and eP = |p|v=Vo and e is a formal small parameter. By considering e 
as a small parameter one may apply to equation (jH]) a statistical perturbation theory that 
was developed by Bogoliubov (1945). Let us briefly formulate the main results related with 
this theory. Let us consider Hamiltonian system with n degrees of freedom, which can be 
described by Hamiltonian H{p,q) in the field of external force ef{t) , where e is a small 
parameter so that complete Hamiltonian can be represented as 

T = H + ef{t)P 

where P{p, q) is some function of canonically conjugated coordinates and momenta. With 
respect to disturbing force we will assume that it is a real function of time like 

fit) = ^ cos{ut + (f^), 

Where phases ipu are independent random variables which are uniformly distributed over 
the circle. Frequency spectrum of disturbing force may be assumed to be almost continuous 
so that the sum will be presented as 

By having continuous Fiu) may be replaced by integrals 

+00 

j F{v)I{v)dv 
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Having dynamical system proposed by Bogoliubov (Bogoliubov, 1945) we may come down 
to description of the main result so that paper will relevant to the case of classic mechanics. 
Let us label probability density for coordinates and momentas via Dt given that phases have 
some defined values (pi,. Probability density for p and q in the usual sense, i.e. under in fixed 
phase values (p^, may be calculated from Df by using averaging operation over all phases. 

At the initial moment of time t — we assume that distribution of coordinates and momenta 
Dq is independent on phase, i.e. 

Do = Po, 

. Evolution of Dt in time should happen according to the well known Liouville equation: 

^{H,Dt) + ef{t){P,Dt), 



dDt 



dt 

where (A, B) is a Poisson Bracket defined as 

,dAdB dAdB 



^ dpi dqi dqi dpi 

Let us introduce another one-parameter group of operators, T^, i e R, that acts on dynamical 
variables according to 

TtF{p,q)^F{pt,qt), 

where 

dH{pt, qt) 



Pt = 
<it = 



dpt 
dH(pt,qt) 



dqt 

with initial data Pq = p, go = <1- In these assumptions and notations it was derived (Bogoli- 
ubov, 1945) (within limits of small s), that the following equation on pt. 



dpt 
dt 



+ 00 

= {H, Pt) J A{t - t){P, Tr^tiP, Pr))dT, (7) 
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where 



A(r) = - I{v)cos{vT)dv. 
2 JO 



It is important to note that as it follows from conclusion of equation ((Tj) (Bogoliubov, 1945) 
this equation might be true if not only phase (p^, but also amplitude for disturbing force 
would have been random values. More correctly, let us assume that disturbing force is 
expressed as More correctly, let us assume that disturbing force is expressed as 



ivt\ 



where, Vz/, u > 0, whereas are random values fiy defined within the same probabilistic space 
{Q, S, P)). Here, Q is some set, S is some algebra of measurable sets, and E is a probability 
(positively defined) measure defined at sets made of so that P{^1) = 1. Let us assume that 
fu satisfies to the following relations 

fu = i: = 0, 



fuffi fufu 



where 



Wfi, u > 0, 
e:= j e(x)rfP(x). 



Frequency spectrum u we will takes practically continuous so that the sums 

J2f{u)I{u) 

V 

with continuous variables F(v) maybe replaced with integrals in the limit 

F{y)l{y)dv 





It is important to note that volume of domain D of cell ^ belongs to a class of random 
variables of such type. Let us assume that at the initial moment of time distribution function 
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depends only on momenta and coordinates of protein molecules in domain but not on 
the point uj of probability space VL. 

Let us again assume pt = Dt ■ Thus, at the initial moment of time t = the following 
will be true pq = Dt. Function pt will continue to satisfy Fokker-Planck equation (5) where 
function A(r) is expressed in terms of /(z/) by using equation ([8]). 

In our case when we speak about a living cell we consider a protein system localized 
inside a domain D, whereas disturbing force is f(t) = V(t) — Vq. 

Now we can mention that Vr G M, T^{P, pt) = (TrP,T^pt) due to the fact that any 
canonical transformation keeps Poisson bracket, i.e. we have 



+ 00 



^ = (H,pt)+e' I Ait-r)iP,{Tr.tP,Tr-tPr))dT, (9) 



Let us show every where in equation (Q that we may consider a value T^P as being inde- 
pendent on time. For this we will represent TtP as Fourier series (integral) 

TtP = Y,bue"'' 

V 

As you may see function VP is additive. The latter means that if we break up domain D 
into subdomains Di ,i = 1,2,3... which intersect each other only on their boundaries, then 

N 

vp = Y,yiP^^ 

1=1 

Where 1^ is a volume of domain Di, and functions Pi are calculated over domains Di in the 
same way as function P is calculated over domain D. Of note, that function Pi depends 
only on momentas and coordinates of proteins in domain Vi . As domain D is macroscopic 
then each of domains Di could be considered macroscopically too, and quantity may be 
chosen arbitrarily large. Let us assume that all domains Di have equal volume V/N so that 

1 ^ 

p = -yp.. 

i=l 

Note that function Pi again may be decomposed into Fourier a series (Fourier integral) 



Then we find Vz/ 



1 ^ 



i=l 
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Let us recall that, domains Di are macroscopic so that interaction between particles from 
different areas may be neglected. Thus, if frequency u is not much close to zero, then 
phases argfo^ for different Di are independent and distributed uniformly over the circle. 
The latter property takes place in situation when frequencies u are not close much to zero, 
therefore = for such frequencies. Thus, value T^-P is practically constant in time, and 
instead of equation ([9]) we will obtain 

+00 

^ = {H,p,) +e' I A{t - r)(P, {P,T,,,p,))dT. 


However, within some small time interval where A(t — r) is significantly different from zero 
we may neglect by interaction in differential equation on pt in the limit of small e (Liouville 

equation), i.e. Tt-rPr = Pt- In other words in limit of small we obtain 

+00 

^ = {H, p,) +e' I A{t - r) (P, (P, p,))dT. 


Finally, in the limit of small e we will have 

^ = (/7,a)+^'^/(0)(P,(P,a)). (10) 

We will denote phase space for proteins in domain D through $ and the set of canonical 
coordinates and momenta for arbitrary point of phase space $ via x. Thus x = {p,q)- 
Of note, in equality f lTOj) we did not put into consideration a summand from complete 
Liouville operator L, that describes dissipation. However, a member describing dissipation 
reflects a systematic drift in phase space and is a first-order differential operator, which 
preserves probability. Thus, final equation for pt, that take into account possible dissipation 
will be like 

^ = Dp, + e'^m{P,{P,p,)), (11) 

where 

2n ^ 

Dptix) = -—ib^ix)pt{x)) 

For some set of real- valued functions &^(x) on phase space of proteins, n is a degree of 
freedom of the system. Now we need to define set of functions b^{x), p = 1,2, ...2n. Let us 
do it by relying on assumption that under equilibrium distribution 

1 h'{x) 

Pt{x) = — e T . 
16 



Right part of equation (fTT]) must become zero. Of note, operator p i— )■ (P, {P, p)) represents 
a second-order real self-adjoint differential operator in L2($, dT), that preserves probability, 
where dT = dpdp is a standard Liouville measure on the space $. Thus, this operator may 
be rewritten as 

2n 

dxi, ^'"^^ ' dx. 



1/ ^-^a 

Where for each x G $ is a real- valued symmetrical matrix a^j,(x). Thus, equation ( ITTj) may 
be rewritten as 

^ =i:ir{^%mf:a,.Mir - (12) 



dt ^ dxu 2 dxy 

(1=1 ^ !/=l 



So we will obtain 



11 = 1 i/ = l 

A particular solution for this equation is 



2n r% 



b,{x) = -e^—I{0)J2a,,{x)—H\x). (13) 



i/=i 



General solution of ( fT3l) is defined modulo the solution of the following equation 

^ -{bUx)e~^}=0 



dxy 

relative to h'^{x). Now let us notice that function /(w), as it will be seen later, may be 
modified by artificial means if change mechanistically an environing ^\D oi area D. Thus, 
a member that describes dissipation in equation ( IT2!) should not depend on /(O). It results 

d,b'^ = 0. 
b'^d.Hix) = 0, 

T ~ /(O) (14) 

We assume that b'^{x) = . However, this assumption will not influence on our further 
conclusions. Finally, there will be 

/i=l u=l 
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where 

^ 2n Q 

b,ix) = -e^l^I{Q)Y^a,,{x)—H\x). 

After introducing preliminary remarks in mathematical terms let us describe how we under- 
stand a mode of action for Ling's medicines. 

Let us introduce two parameters Ax and A2 , where parameter Ai represents concentra- 
tion of medicine inside cell, and parameter A2 is a concentration of protein molecules that 
bound molecules of the medicine. In physical terms it is clear that the increasing of Ai leads 
to increasing of A2 (until molar concentration of the medicine will be lower compared to 
molar concentration of the target protein), so we will focus on A2 only. 

Then, let us assume that concentration of medicine inside cell is going up, which is 
accompanied with the corresponding increase of A2. Binding of one molecule of medicine 
to a protein molecule leads to redistribution of electron density over protein (Ling, 2001). 
According to Ling, it makes protein molecule more hydrophilic (due to destruction of the 
secondary structures) that adsorbs on its surface more water. Formation of quasi- crystalline 
adsorption layer provides protein/water complex with elastic properties. 

Let us consider a thin fiat homogenous metal plate (e.g., made of iron) that has temper- 
ature T, and small, but macroscopic area D on it. We will examine fluctuations of volume 
D. Assume that spectral density for intensity of volume D fluctuations is defined as /(w), 
a; G M, > 0. Now, let us assume that we marked by disk an area on this plate A that does 
not intersect with D, and then did the following: we cut off area of the plate and welded 
instead of it another disk of the same size, but made of more elastic metal, e.g. titanium. It 
seems obvious that after such manipulations spectral density of intensity for volume should 
be reduced. In our case a role of "plate" will be played by an injury area of the cell. If 
protein of such area start to bind medicine then it will lead to formation of physioatoms at 
resting state which may be identified as "disks" S (physioatoms and their associates), which 
have higher elasticity ingrained into an area of injury. We need to mention that equilibrium 
distribution when the right part of Fokker-Planck equation (|T2|) is being set to zero will be 
as 

I h'(x) 

Where effective temperature Tgjj is linked to spectral density of intensity at zero /(O) by 
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the following proportion 

In other words 

Teff = const/ (o) 

Thus, finally binding of medicine to a particular protein molecule results in reduction of 
effective temperature T^ff for proteins in area D. It means that one molecule of medicine 
may exert therapeutic action on several protein molecules. 

While discuss in generalized thermodynamics in Part 2 we talked about existing critical 
value for temperature Tcr, above which no life can exist. In ability to maintain life above 
temperatures Tcr we associated with the fact that entropy looses plateau as it is a function 
of values of non-trivial first integrals commuting with each other under fixed energy level. 
The consequence of this is that difference in description between generalized microcanonical 
distribution is used (that we believe describes living entity at resting state) and conventional 
equilibrium distribution will disappear. In other words difference between animate and 
inanimate will disappear. 

When proteins adsorb extra amount of water while being acted upon by medicine then it 
means that a medicine hinders processes of excitation/injury, and facilitates transition of the 
cell towards resting state. Acting together with ATP a medicine contributes to destruction 
of non-native protein structures, protein aggregates dissociate, individual proteins unfold, 
and the consequences of damage disappear. Medicine and serves as an independent factor, 
and as synergist of ATP. 

Additionally, we can make another point, although hit does not have direct connection 
to theory about mode of action of medicines, but rather represents a self-contained interest. 
We assumed that below critical temperature state of protein system may be described by 
generalized microcanonical distribution (see Part 2), and if to describe system by using such 
distribution will not be the same as if it was described by standard equilibrium distribution. 
At the same time, we found that stationary solution for Fokker-Planck equation f ll2p gives 
a regular Gibbs equilibrium distribution. However, we believe that our generalized micro- 
canonical distributions may be obtained by applying Bogoliubov's quasi-average approach. 



2n 



—e 



2T 



d 



/(0)5^a,.(x)— i/'(x). 



i/=i 
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For example, recently (Prokhorenko, 2011) there was shown that, on one hand, supercon- 
ductive state in Bardeen-Cooper-Schrieffer theory of superconductivity may be reached by 
using Bogoliubov's quasi-average approach, but, on the other hand, by applying generalized 
thermodynamics. Roughly speaking, we think that when system is described by generalized 
microcanonical distribution giving plateau for entropy (as a function of values for non-trivial 
commuting first integrals under certain energy) then it means that it is not merely similar 
with the same condition in theory of phase transitions, but truly represents a typical case of 
multi-phase area. Also, we believe that our generalized thermodynamics may be restated in 
terms of theory of phase transitions if the latter will be properly developed. 

4 Special features of ATP as Ling's medicine 

Despite the fact that from the very beginning Ling's medicines have been defined as ATP 
analogues based on their mode of action (electron-accepting compounds), however, in context 
of the theory that wc propose, properties of ATP are not considered so unequivocal compared 
to medicines which are administered from the outside. 

In the previous part we proposed two parameters Ai and A2 for concentration of molecules 
of medicine inside cell and for concentration of protein molecules which are controlled by 
medicines, respectively. In the previous part we proposed two parameters Ai and A2 for 
concentration of molecules of medicine inside cell and for concentration of protein molecules 
which are controlled by medicines, respectively. These two parameters may be applied for 
ATP as well provided that parameter A2 is an increasing function for parameter A2: when 
Ai grows up so does parameter A2, i.e. a number of proteins which are controlled by ATP 
increases too. Variability of states (fluctuations) for physioatoms and their associates goes 
down due to increased stiffness of target structures for ATP. Thus, effective temperature of 
protein molecules in the cell decreases whereas cell gets rid of damaged structures. 

Here, however, there is another question. In the theory of mode of action for medicines 
(Part 3) we consider Ai (which means A2 as well) as a parameter, that may be changed in 
an arbitrary way. However, in case of ATP we can not further change Ai at our wish, as 
ATP is being synthesized inside the cell and depends on the cell functional condition. It 
rises a question: whether our theory about mode of action for medicines may be applied for 
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ATP? 

Contradictions may appear on several particular points. 

1. On one hand, we showed that during Ling's cell death its key proteins get folded 
(Prokhorenko and Matveev, 2011a,b). On the other hand, our calculations claim that de- 
creasing of ATP concentration Ai under excitation/injury results in increase of effective 
protein temperature. Because increase of temperature enhances the probability that protein 
will exist in unfolded conformation, then decrease of ATP concentration seemingly have to 
result in increased number of unfolded proteins inside the cell. 

However, its hold be kept in mind that effective protein temperature is not equal to 
temperature of cellular water, meaning that when effective protein temperature grows up, 
so when ATP concentration decreases, then temperature of cellular water will be kept un- 
changed. Generally, under increase of cell temperature proteins start to unfold, mainly due 
to increased intensity of hits on protein molecules performed by water molecules. However, 
as in the context of our analysis, temperature of water does not change, it means that it 
does not contribute to increased half-life for unfolded state. 

By using Le Chatelier principle (Landau and Lifshitz, 1995) we may define conditions 
(factors) which may cause folding of "hot" protein molecules. Let us consider this issue in 
more detail. 

Generally. Le Chatelier principle is defined as: external force that disturbs equilibrium 
of an object stimulates processes inside it, which counteracts the force to weaken results of 
its influence. However, we will need to consider only special case of this principle, so now 
we will get down to the correct definition of it. 

Let us consider some entity under constant pressure and let y he a parameter that 
describes it and characterizes when an entity is at equilibrium. As a system is considered 
under constant pressure then a condition of its equilibrium may be written as 

^ d^P,T,y) ^ ^ 
dy 

Where $(P, T, y) is a thermodynamic potential. Then, it will be correct to write an inequahty 
(Landau and Lifshitz, 1995). 



,dT dT 



y > i^)y > (15) 



Notably, in the latter formula both derivatives were calculated under constant pressure. In 
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our special case for living cell as a parameter y we will use effective temperature of protein 
molecules Te//. In this case condition F = may be equivalently expressed by equation 
Teff = T. From equation (fTSl) there will be: 



Thus, we may conclude 



(^)t > 0. (16) 

Let parameter characterizing a degree of protein unfolding, that was previously 

introduced (Prokhorenko and Matveev, 2011a) so that proteins are considered to be more 
unfolded when gets bigger x. Then, from equation (flGll there will be obtained: 

However, if under constant pressure and temperature proteins will start to unfold under 
influence of any external action, then they will start to adsorb water so that according to 
Ling entropy of water will be decreasing. The latter statement is relevant and true as while 
water gets absorbed and the number of microscopic states which represents this particular 
macroscopic state, will be decreased. In other words, by Ling 

Thus, from inequality f fTTl) there follows that 

, dx , 

n < 



meaning that when effective temperature of proteins Te// increases under unchanged tem- 
perature T of cellular water then protein molecules get folded. This is what we wanted to 
find. 

2. When cellular ATP concentration goes up Ai, it means that a number of physioatoms 
in the cell will be increase in gas well, including their size (due to more effective adsorption of 
water) and a number of associates of physioatoms. When an amount of quasi-crystal water 
inside the cell (stabilized by unfolded proteins) goes up, then stiffness of the cell will grow 
as well. It rises a question: what is possible to say about changes of module of Ling's cell 
isothermal compressibility (under constant pressure) during cell death (when its stiffness 
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decreases) basing only on generalized thermodynamics? If it turns out, that according 
to generalized thermodynamics module of isothermal compressibility grows up instead of 
decreasing then it will contradict to corollaries from our theory of Ling's medicines. 

Now we will show that generalized thermodynamics results in statement that under cell 
death module of isothermal compressibility goes down. Thus, let introduced in Part 2 will 
be a continuous parameter that describes a number of active involutive first integrals in cell. 
It is relevant to remind that a condition corresponds to situation when all first integrals 
of interest are active, whereas condition corresponds to situation when none of integrals of 
interest is active. 

Thus, let s G [0, 1] introduced in part 2 will be a continuous parameter that describes a 
number of active involutive first integrals in cell. It is relevant to remind that a condition 
s — corresponds to situation when all first integrals of interest are active, whereas condition 
s = 1 corresponds to situation when none of integrals of interest is active. 

Let s ^ s — Ss , where (5s > is an infinitely small quantity. Let us assume that 
f{S) = (AS) E,T is a change of cellular entropy happening during this process under constant 
energy and volume. Now we will have 

f{S) = {AS)e,v = {AS)t,v - %i^S)T,v = {AS)t,v - ^{AS)t,v (18) 



Further 



that results in 



^/(S) = ^((A5)Ty - i(AS)T,v) = ^(AB)Ty (19) 



(^) = ^(AE)ry (20) 



•^/(^)) . (i),(^), (21) 



dS ' ^dS'"^ dT 



Next, there will be 



^^^^{AE)tv. (22) 

There exists such a temperature Tq, that for T > Tq {AE)Ty = animate state of matter 
will not be possible to distinguish from dead matter, and under {AE)Ty — 0. Thus, under 
T < To, {AE)Ty < if value of entropy S corresponds to value of temperature T. 
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Let us consider operator (^)t that acts only on function of entropy of system. Then, 



>9P 

we will have 



, 9 , .dS , . d , /^^\ /^^ \ ( ^ \ (no\ 

i-Q^h = ig^M-Qsh = (dpM-g^M-Qsh- (23) 

But the following thermodynamic inequality states that 

dV 

(^)t < 0. (24) 

However, the majority of entities under adiabatic expansion will get cooled so that we may 
assume that 



Thus, 



(99 

igyh > 0. (25) 



igphigyh < 0. (26) 

and 

Qt < (27) 

Further, parameter s characterizes a number of first integrals in involution. Let he s ^ s—Ss, 
where SS > is an infinitely small quantity. We may remind, that f{S) — {AS)e,t- 
Then under certain energy and due to the fact that (|f )y — T, there will be stated that 
under s^s~5s^E^E — Tf{S). Thus, according to theorem about small additives 
F 1-^ F — Tf{S), where F is a free energy, and $ i-^ ^ — Tf[S), where $ is a thermodynamic 
potential. According to the following formula 

(9$ 

V = igph. (28) 



Module of isothermal compressibility 



dV 9^$ 

kr := {qp)t = (^)t (29) 



If s !->■ s — (5s , then kx — kx — Skx, where 



^2$ 

SkT = T(^)^ (30) 
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Further, we will have, 



But 



As a result 



r((|.).(|.)./(S))^T((A),(|^),fm) (31) 



dS (9^$ d (9$ dV 

^dp'^^ ^ ~dPdT ^ ~&rdP ^ "^ar^^ ^^^^ 



(33) 



In other words 



- -n-Q^fis) + {^)p{^)Tns)} (34) 

After that we will find that 

= + (35) 

Eventually, we find that the resultant expression for change of module of isothermal 
compressibility is 

rfV r)V (99 

Further, -^f{S) = j^{AE)t,v and f'{S) — ii S > Sq for some 5*0. A cell can not 
lives under temperatures above certain critical value and under temperatures below certain 
critical value {AE)t,v < because there was proved (Prokhorenko and Matveev, 2011a) 
that the cell when it dies under constant volume and temperature it releases energy, so 
death is an exothermic process. Thus, if f'{S) is not too "pathological" function, then there 
could exist such 5" < 5*0 that f'{S) increase within interval [S",S'o]. So, within interval 
[5", So] f"{S) > 0. Further, there exists such value 5"' < So , that function f'{S) in the first 
summand may be neglected as j^f{S) = j^{AE)t,v = if ,5 > S". Thus 

^ -Q;)p{^)Tf"{S) (37) 

However, we decided that (|p)t < 0. Further, we may state that (§^)p as the majority 
of entities when they are heated they get expanded. Therefore, under these assumptions 
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Skx > 0) , but under s ^ s — Ss ^ — Skx, i.e. absolute value of module of isothermal 
compression under constant pressure will increase and the cell will become less stiff. This is 
what we wanted to prove. Likewise, one may say that under cell death absolute module of 
isothermal compression at constant pressure will increase as well. 

5 Conclusion 

In the current paper we discussed a probable physical mechanism that describes a mode of 
action for medicines belonging to electron-accepting compounds. Basically, our approach 
relies on so-called generalized thermodynamics (Prokhorenko and Matveev, 2001 a). An ex- 
plicit physical interpretation of this mechanism is as follows. When concentration of medicine 
inside the cell goes up it results in increased number of protein-medicine complexes. After 
that electron densities on key functional groups of protein molecules are being redistributed 
that leads to adsorption of water by target proteins and further expansion of adsorbed water 
layer. Reduced motility of water molecules in adsorbed layer results in situation when effec- 
tive temperature of protein molecules drops down as well all processes in protoplasm which 
are connected to excitation and damage stop, i.e. the cell gets recovered. 

Further, we used the proposed theory for mode of action of medicines that may be applied 
to explain a role of ATP in cell vital activity in context of Ling's theory. We consider that 
there is a principal similarity of a certain kind of compounds (Ling's medicines) and ATP. 

We are very grateful to P. Agutter, A.V. Koshelkin, Yu. E. Lozovik, A.V. Zayakin, E.N. 
Telshevskiy for valuable critical comments on this article and very useful discussions. 
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